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Abstract 

We study the stability problem for a non-relativistic quantum system in dimension three com- 
posed by N > 2 identical fermions, with unit mass, interacting with a different particle, with 
mass m, via a zero-range interaction of strength a G M. We construct the corresponding renor- 
malised quadratic (or energy) form T a and the so-called Skornyakov-Ter-Martirosyan symmetric 
extension H a , which is the natural candidate as Hamiltonian of the system. We find a value of 
the mass m*(N) such that for m > m*(N) the form J- a is closed and bounded from below. As 
a consequence, T a defines a unique self-adjoint and bounded from below extension of H a and 
therefore the system is stable. On the other hand, we also show that the form T a is unbounded 
from below for m < m*(2). In analogy with the well-known bosonic case, this suggests that the 
system is unstable for m < m*(2) and the so-called Thomas effect occurs. 



1 



1 Introduction 



The dynamics of a quantum system composed by N particles in R d , d = 1,2,3, interacting via a 
zero-range, two-body interaction is described by the formal Hamiltonian 

N . N 

n = "E^- A ** + E tMjS(*i-Xj), (1-1) 
1=1 * ij=l 

«i 

where Xj G R rf , i = 1, . . . , N, is the coordinate of the i-th. particle, mi is the corresponding mass, A x . 
is the Laplacian relative to Xj, and G R is the strength of the interaction between particles % and 
j. To simplify the notation we set h = 1. Formal Hamiltonians of the type (jl.ip are widely used in 
physical applications. In particular they are relevant in the study of ultra-cold quantum gases, both 
in the bosonic and in the fermionic case, in the so-called unitary limit, i.e., for infinite two-body 
scattering length (see [BH| . [CWl] . [CW2| . [CMP] and references therein). 

The first step towards a rigorous approach to the analysis of the model is to give the mathematical 
definition of such a Hamiltonian as a self-adjoint operator on the appropriate L 2 -space. One first 
notices that the interaction term in (jl.ip is effective only on the hyperplanes Uj<j{xj = x^}. This 
suggests to consider the operator %q defined as the free Hamiltonian restricted to a domain of smooth 
functions vanishing in the neighbourhood of each hyperplane {xj = xj}. It is easily seen that Ho 
is symmetric but not self-adjoint and a trivial self-adjoint extension is the free Hamiltonian on its 
natural domain. Then, by definition, any non trivial self-adjoint extension of the operator %q is a 
Hamiltonian for a system of N quantum particles in R rf with a two-body, zero-range interaction. 

The second and more relevant problem is the concrete construction of such self-adjoint exten- 
sions. It turns out that each extension is characterized by a specific generalized boundary condition 
satisfied by the wave function at the hyperplanes. The two most frequently used techniques for 
the construction are Krein's theory of self-adjoint extensions and the approximation by regularized 
Hamiltonians, in the sense of the limit of the resolvent or of the quadratic form. The difficulty of 
the analysis depends on the dimension d. For d = 1 the problem is greatly simplified by the fact 
that the interaction term is a small perturbation of the free Hamiltonian in the sense of quadratic 
forms. For d = 2 a natural class of Hamiltonians with local zero-range interactions was constructed 
in [DFT] exploiting renormalised quadratic forms and it was also shown that such Hamiltonians are 
all bounded from below. For d = 3 the problem is more delicate. In order to illustrate the main point, 
let us consider the special case n = 2 where, in the center of mass reference frame, one is reduced 
to study a one-body problem in the relative coordinate x with a fixed <5-interaction placed at the 
origin. In this case the problem is completely understood (see, e.g., [AGH-KHJ) and the entire class 
of Hamiltonians can be explicitly constructed. It turns out that the domain of each Hamiltonian 
consists of functions ip G L 2 (R 3 ) n if 2 (R 3 \ {0}) which exhibit the following singular behaviour for 
|x| -»■ 

q 

ip(x) = — r + r + o(l), with r = aq, (1.2) 

where q G C and a G R is a parameter proportional to the inverse of the scattering length. We un- 
derline that the relation r = aq in f)1.2|) should be understood as the generalized boundary condition 
satisfied at the origin by all the elements of the domain. In the general case N > 2 the characteri- 
zation of all possible self-adjoint extensions of Hq is more involved. However, a class of extensions 
based on the analogy with the case N = 2 can be explicitly constructed. More precisely, one con- 
siders the so-called Skornyakov-Ter-Martirosyan (STM) extension of Ho which, roughly speaking, is 
a symmetric operator acting on functions ip G L 2 (R 3Ar ) n H 2 (M? N \ Uj<j{xj = x,}) satisfying the 
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following condition for |xj — Xj| — > 0: 

ip(x.l, ■ ■ ■ ,x n ) = | — ^ — - + + o(l) , with nj = aijqij , (1.3) 

where qij is a suitable function defined on the hyperplane {x, = xj} and {a^} is a collection 
of real parameters labelling the extension. Noticeably, the boundary condition (jl.3p defining the 
STM extension of Ho is a natural generalization to the case N > 2 of the condition (|1.2|) that 
characterizes the two-body case. Unfortunately, unlike (jl.2p . (jl.3p does not necessarily define a self- 
adjoint operator. Indeed, for a system of three identical bosons it was shown in [FM] that the STM 
extension is not self-adjoint and all its self-adjoint extensions are unbounded from below owing to 
the presence of an infinite sequence of energy levels going to —00 for k — > 00. In [MM] this 
result was generalized to the case of three distinguishable particles with different masses. This kind 
of instability is known in the literature as the Thomas effect. It should be stressed that the Thomas 
effect is strongly related to the well-known Efimov effect (see, e.g., [BHj . [ADFGL] . |AH-KW] ) even 
if, to our knowledge, a rigorous mathematical investigation of this connection is still lacking. We also 
mention that if, instead of (|1.3p . one introduces a "non- local" boundary condition on the hyperplanes 
then it is possible to construct a positive Hamiltonian and to study its stability properties for N 
large (see, e.g., jFSI). In this paper we do not consider this kind of Hamiltonians. 

It is reasonable to expect that the Thomas effect does not occur if the Hilbert space of states 
is suitably restricted, e.g., introducing symmetry constraints on the wave function. A remarkably 
important constraint is antisymmetry. In fact, a wave function that is antisymmetric under exchange 
of coordinates of two particles necessarily vanishes at the coincidence points of such two particles, 
thus making their mutual zero-range interaction ineffective. Analogously, in a mixture of fermions 
of different species subject to pairwise zero-range interaction, fermions of the same species cannot 
"feel" mutual zero-range interaction and therefore the interaction term in the Hamiltonian is less 
singular. 

In this paper we consider the simplified model consisting of N identical fermions, with unit mass, 
and a different particle with mass m, interacting with the fermions through a zero-range potential. 
For such a model only partial results are available and it is remarkable that they strongly depend on 
the parameters N and m. 

Concerning the physical literature, we mention that for N = 2 it is known (see, e.g., |BH| 
and references therein) that for m < 0.0735 = (13.607) -1 the Thomas effect is present while for 
m > 0.0735 the STM extensions are expected to be bounded from below. More recently (see [CMP]), 
it was shown by means of analytical and numerical arguments that in the case N = 3 the Thomas 
effect occurs for m < 0.0747 = (13.384)" 1 . This, in particular, indicates that for 0.0735 < m < 0.0747 
there is a sequence of genuine four-body bound states with energy going to —00. 

From the rigorous point of view, the case N = 2 was investigated first in [Ml] and |MS| . where 
it was proved that the STM extension is self-adjoint if m = 1. In [S], the existence of a critical 
mass m*(2) ~ 0.0735 was shown such that for m < m*(2) any STM extension - more precisely 
its restriction to the subspace of angular momentum I = 1 — is not self-adjoint and its self-adjoint 
extensions are unbounded from below. Therefore, the system is unstable and the Thomas effect 
occurs. This result was extended in [FT], where it was shown that the quadratic form associated 
with the STM extension, restricted to the subspace of angular momentum I, is unbounded from 
below if an explicit condition on m and I is satisfied. Finally, in the case A^ < 4 and m sufficiently 
large it was shown in [M3| that the STM extension is self-adjoint. 

In the present work we study the problem for generic A" and m and following the line of |DFT] 
we construct a renormalised quadratic form T a , a£l, which is naturally associated with the STM 
symmetric extension H a . Here a is the value that all the parameters aij labelling the STM extension 
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must be equal to as a consequence of the fermionic symmetry. Note that in the bosonic case the 
quadratic form would differ in a sign in front the non-diagonal term (see the remark after (|2.15p ). 

The first question we address is a non-trivial sufficient condition for the stability of the model. 
As a first main result (Theorem 12 .2p we prove that for any N there is a value of the mass m*(N) > 
such that J- a is closed and bounded from below if m > m*(N). This implies that T a is the quadratic 
form of a unique self-adjoint extension of the STM operator H a , and this extension is bounded from 
below. It therefore describes a stable system, where the Thomas effect does not occur. 

Such a critical mass was first conjectured in [M2] and is precisely the unique root of an explicit 
equation (see (|2.16p and (|2.18j) in Section [2]). It turns out that m*(N) is increasing with N and 
that the condition m > m*(N) guarantees the stability also in the limit of infinitely many fermions, 
provided that the mass of the extra particle scales asmociV. 

The second question we address is a sufficient condition for the instability of the model. This can 
be seen by plugging suitable trial functions into J- a . An attempt in this direction is in [DFT1 Section 
7], but with trial functions that do not satisfy the fermionic symmetry: thus, the result stated there 
on the unboundedness from below of the quadratic form for m = 1 and N sufficiently large cannot 
be considered valid. Our second main result (Theorem 12. 3D fills in this gap and we prove that for 
any N > 2 the quadratic form T a is unbounded from below for m < m*(2). In analogy with the 
bosonic case, we expect that in such a case H a is not self-adjoint and all its self-adjoint extensions 
are unbounded from below. 

Let us make a few remarks on the above-mentioned results. First, we emphasize that in the case 
A = 2 we fully characterize stability: the model is stable for m > m*(2) and unstable m < m*(2). 
Whereas the latter was already found in [Sj by means of the theory of self-adjoint extensions, the 
former is proved here for the first time. 

Second, if N > 2 we expect the condition m > m*(N) to be far from optimal for stability. This is 
due to the crucial role played by the restriction to antisymmetric wave functions (see the discussion 
in Section U|) so that the system might be stable also if our condition is violated. 

Third, the fact that the (A + l)-particle system is unstable at least when m is below the same 
threshold m*(2) for the instability of the (2+l)-particle system has a rather natural interpretation: 
the instability of a subsystem made of two fermions plus the different particle is responsible for the 
instability of the whole system. 

We want to mention that in the final stage of the preparation of this work we became aware of a 
recent paper [M4j where the case of two fermions plus a different particle is studied using the theory 
of self-adjoint extensions. We believe that a comparison with our methods and results would be of 
great interest for the further developments of the subject. 

The paper is organized as follows. In Section [2] we introduce the renormalised quadratic form J- a 
and the STM extension H a and we formulate our main results. In Section [3] we give the proof of 
Theorem 12.21 In Section U] we give the proof of Theorem 12.31 In the Appendix we briefly outline the 
formal renormalisation procedure to derive J- a . 

For the convenience of the reader, we collect here some useful notation that will be used through- 
out the paper. We use the notation Lf(R. d ) (resp. H^(M. d ), H { 1 ^ 2 (M. d ), etc.) for the space containing 
totally antisymmetric functions belonging to L 2 (M d ) (resp. 1 (IR. d ) , H~ 1 ^ 2 (M d ), etc.). We often use 
the short-hand notation || • \\ L 2, \\ ■ \\ L 2, etc. for the associated norms || • H^^d), || • ||x, a (]R d )5 e * c - 
For a vector x £ I 3 we set x = |x|. Moreover we define K := (k.%, . . . , kjy-i) and, for i = 1, . . . , N, 

Kj := (ki, . . . , kj_x,kj + i, . . . , kjv) • 

For any / G L 2 (R d ) the Fourier transform is defined by /(k) = (2-7r)~ o! / 2 J Rd dx e _lk ' x /(x) . 
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The functions G x , Gx, L\, with A > 0, and D(K) are defined in p?7j) . (1231) . (pT8l) and fl33J) respec- 
tively. 



2 Main results 

In this section we introduce the quadratic form J- a , the STM extension H a and we formulate our 
main results. 

2.1 The quadratic form T a 

For a three dimensional quantum system composed by N identical fermions, with mass one, plus 
a different particles, with mass to, with a two-body zero-range interaction the formal many-body 
Hamiltonian is 

1 1 N N 

^ := " 2ro" Axo " 2 ^ Axi + ^ ^ 5(X ° " Xi) ' (2-1) 
i=i i=i 

where Xj 6 M 3 , i = 0, . . . , N, and \i £ M. Introducing the centre of mass and relative coordinates 

1 ^ 

X := — ( tox + ^ x 

to + N V ^ 

i=l 

k yi:=x -Xj, for i = l,...N 



m + N y—v 1 z^v (2.2) 



one obtains 



H = H cm + ?g±H (2.3) 



cm ^ 2m 

-1 



here H cra := -[2 (to + iV)]" 1 Ax and 

^ : = - E A y> " ^TT E v y* • v ^ + m E , (2-4) 

i=l i<jf i=l 

V y . denoting the gradient with respect to yj. We also introduce the free Hamiltonian 

N 2 

S(ff ) := HhR 3N ), H := - £ A yi - _ £ V y , • V y , , (2.5) 

i=l i<j 

and its restriction to functions vanishing in a neighbourhood of the hyperplanes {yi = 0} 

9{Hq) :=Ul) e #f (R 3iV ) jfdk i ^(ki J ...,k JV ) = 0, i = l,...,iv|, #„ := #o| s( £ o) • (2.6) 

In order to give a rigorous meaning to the formal expression (12, 4h as a self-adjoint operator in 
Lj(lR 3Ar ), one can use Krein's theory to construct self-adjoint extensions of the operator (|2.6p (see 
|M3| ). Here, instead, we follow a different approach and we investigate the quadratic form associated 
with the expectation value (ip\H\ip). However, because of the singularity of the "potential" <5(yi) 
the quadratic form has to be defined via a renormalisation procedure in the Fourier space. The idea 
of the construction is given in the Appendix (see also [DFT],[FT|) and here we only give the final 
result. 
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Let us denote for any A > 

G A (ki, . . . ,k N ) 



N 



Ekf H k { ■ k, + A 

i=l i<j 



(2.7) 



L A (k l5 . . . , kw _0 := 2vr 2 f E ^ 2 + -^- 2 k, • k, + A) , (2.8) 

V (m+iy ^ (m + 1) 2 J 

and, for any "charge" £ G .ff f 1//2 (M 3iV ~ 3 ), let us denote the "potential" produced by £ in the Fourier 
space by 



AT 

&C) (k x , . . . , kjv) := ^(-l) l+1 G A (k!, ...,k N ) i(Ki 



(2.9) 



»=i 



where Kj := (kx, . . . , k_i, k+ij • • • ; kv). Then the quadratic form J 7 ^ in the Hilbert space L 2 (IR 3Ar ) 



is defined as follows 



f 2 (M 3JV ) 



3£g^(<&*) s.t. 4 -g^sH^R 3 ")} 



2 mR3N) + iv$^ [e] , 



where A > 0, J-} 



0[<pj := (0| -ffo |0) an d ^ is the following form on the charge £ 

: =i/ f 1/2 (M 3Ar - 3 ), 

*aKl :=< a A S [e]+*f[C], 

dki • • • dkjv-x £(ki, . . . ,kv-i) [a + L A (ki, . . . , k/v-i)] 



(2.10) 
(2.11) 

(2.12) 
(2.13) 
(2.14) 



$f [£] := / dsdtdk 2 ---dk JV _ 1 |*( 8 ,k 2 ,...,k JV _i)|(t,k 2 ,...,k JV _i). 

JM? N 

G A (s,t,k 2 , . . . ,kjv-i^ 



(2.15) 



We notice that for £ G @(&a) we nave £a£ G Lj(R 3N ) and £ A £ ^ ^(IR 3 ^). Therefore the decompo- 
sition in (|2,10p is meaningful. Moreover, as we shall see in Section [3l the form <I>^ is well-defined on 

no- 

It is worth mentioning that the fermionic constraint implies not only that the wave function is 
totally antisymmetric but also that the form on the charges f)2 . 13|) differs from the bosonic case by a 
sign in front of the off-diagonal part (|2.15p . This fact results in a weaker effective interaction among 
the fermions and the stability problem is qualitatively different. 

In order to formulate our main results on the form T a we first introduce our definition of stability 
parameter. Let us consider the following function 



A(m,N) :=2tt~ 1 (N - l)(m + if 



1 



^m(m + 2) 



arcsm 



1 



m + 1 



It is easy to check that, for N fixed, the function A(m, N) is decreasing in m and 

lim A(m,N) = oo, lim A(m, N) = . 

m— s>0 m— >oo 

Then we have 



(2.16) 



(2.17) 



6 



Definition 2.1 (Stability parameter m*(N)). 

For N fixed, we define m*{N) as the unique solution to the equation 



A(m,N) = l. (2.18) 

If we define 6 := arctan yj\ + ^, where 6 G (f > §), a direct computation shows that equation (|2.18p 
can be equivalently written as 

cot 26 + 26 - | ^1 - j^-j cos 2 29^j = . (2.19) 

We remark that (|2.19p reduces for N = 2 to the equation found for the critical mass in |PQ p. 
12871, note 10]. We also notice that m*(N) is positive and increasing with N. In particular, the 
condition A(m,N) < 1 is equivalent to m > m*(N). This condition is crucial to guarantee closure 
and boundedness from below of the form T a . 

Theorem 2.2 (Stability for m > m*(N)). 

Let N > 2 and m > m*(N). Then the quadratic form T a is closed and bounded from below. In 
particular, it is positive for a > and 

for a < 0. 

The proof will be given in Section [3] 

Concerning the instability problem, our result is the following. 

Theorem 2.3 (Instability for m < m*(2)). 

Let N > 2 and m < m*(2). Then the quadratic form is unbounded from below for any a£R, 

This theorem generalizes the result obtained for N = 2 in |FT] and it is proved in Section [H In the 
case N = 2, Theorems 12.21 and 12.31 show that the system is stable for m > m*(2) and unstable for 
m < m*(2), therefore our analysis is complete. When N > 2 the problem remains open: no rigorous 
result is available for m*(2) < m < m*(N). 



2.2 The STM extension H a 

Here we recall the standard definition of the STM extension H a and its connection with the quadratic 
form T a . Moreover, as a direct consequence of Theorem 12.21 we prove our main result on H a . 

As we mentioned in the introduction, H a is a distinguished symmetric extension of the operator 
(|2.6p . For the details of the construction we refer to [M3] and here we only give the definition 

3£ G H* /2 (R 3N - 3 ) s.t. (j>x:=^- Gxt € H?(R 3N ), 

[ dk i ^ x {k 1 ,...,k N ) = (A^ i i)(K i ),i = l,...,N), (2.21) 
(H a + X)ip = (H + A)0 A , (2.22) 



@{H a ) := <j V G LUM. 3N ) 
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where A > and 

JV 



a + LxCKn 



i(Ki)- v (-iy +1 [ dk i c? A (k 1 ,...,k jV )e(K i ). (2. 



23) 



The last equality in (|2.2ip should be understood as the boundary condition satisfied by any ifi £ 
3>{H a ). In fact, by a straightforward computation, one verifies that such equality implies the follow- 
ing asymptotic condition for R — > oo 

[ dk^(ki,...,kjv) = 47r J R(-l) i+1 e(Ki) + a(-l) m |(K i ) + o(l) (2.24) 

Jki<R 

and this is exactly the standard boundary condition, formulated in the Fourier space, characterizing 
the STM extension (see [M3] ). 

The next step is to establish the connection with the form T a . For any ip € @(H a ) one has 

(V| H a + A |V> = (VI H + \ \<f> x ) = (<f> x \ H + A |0 A ) + (Gxt\ H + X \<f> x ) 
= ( ( f ) x\H + X\4> x ) + V(-l) i+1 / dKiCiKi) f dk^ A (k 1 ,...,k w ) 

~{ 7R3JV-3 J R 3 

N 

= {fa] H + X \cf>x) + V(-l)* +1 / dK, tiKi^A^MKi) . (2.25) 
Exploiting (|2.23|) and the antisymmetry of £ we have 

TV 

j>i) i+1 / dK,r(K,)(^i)(K 4 ) 

= J2 akiittkitfia+Lxifti)]- V(-ir +i / dk 1 ---dk JV e*(K i )e(K,)G A (k 1 ,...,k J v) 



= N*?g[Z]+N{N-l)*f[t]. (2.26) 

Therefore we conclude 

(VI H a |V) = T a [V] , V e . (2.27) 

From the above relation and Theorems I2.2| 12. 3| we obtain the following result. 

Theorem 2.4 (Self-adjointness and boundedness from below). 
One has the following two possible alternatives: 

(i) If N > 2 and m > m*(N), then T a defines a unique self-adjoint and bounded from below 
extension H a of the operator H a . In particular, H a is positive for a > and 

mf ^> a - Ml-°AKA0y f0ra< °- <2 - 28) 

(ii) If N > 2 and m < m*(2), then, for any a £ K, no self-adjoint extension H a of H a such that 
S>(H a ) C &{J- a ) can be both self-adjoint and bounded from below. 
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Proof. The assertion (i) immediately follows from Theorem 12.21 and the variational characterization 
of the infimum of the spectrum of a self-adjoint operator. 

The converse (ii) can be obtained by exploiting the Birman-Krein-Vishik theory of positive self- 
adjoint extensions (see, e.g., [ASi], |P1| . [P2] ). We omit the details and refer to \FT\ Proposition 4.1] 
where an analogous results was proven. □ 



We remark that, in the case N > 2 and m < m*(2), an analogy with the bosonic case suggests that 
any self-adjoint extension of H a is unbounded from below and the Thomas effect occurs, although 
such a complete characterization is beyond our purposes. The most natural tool to approach the 
analysis of this case is indeed the theory of self-adjoint extensions (see, e.g., [S|,[M4]). 
We conclude this section with a brief comment relative to the special case N = 2, m > m*(2). In 
particular we want to give the explicit characterization of our Hamiltonian H a in order to make more 
transparent a possible comparison of our result with the result in 



In Section 3 we shall prove the estimate (I3,52h . which in particular implies that the form defined 
on i? 1/,2 (lR 3 ), is closed and positive for A sufficiently large. Therefore it defines a positive, self-adjoint 
operator r*, @{T X ) in L 2 (M 3 ) explicitly given by 



#(T£) ={e E # 1/2 (M 3 ) I T x a i £ L 2 (M 3 )} , 
(r x j) (q)= [a + L A (q)]£(q)+/ dp G A (p, q)|(p). (2.29) 
Exploiting this fact and following the same line of [DFTl Section5], we obtain 



3t£®(T x a ) s.t. 4> x := V-<?a£e tf f 2 ( 



/ 



d P 0A(p,q) = (r^)(q) , (2.30) 



(H a + \)ifi = (H + \)<j> X . (2.31) 

We notice that the above operator differs from the operator (I2.21h . (|2.22p only in a larger class of 
admissible charges £, i.e., the domain @(T X ) strictly contains H 3 / 2 (R 3 ). We also underline that 
the boundary condition satisfied on the hyperplanes by an element of (|2.30p is the standard STM 
boundary condition. 



Closure and boundedness from below of T, 



a 



The proof of Theorem 12.21 is based on a careful estimate from below and from above of the form 3> x 
on the charge £. If > 2 (with an obvious modification in the case = 2) we rewrite both the 
diagonal and the off-diagonal terms of defined in (|2.14p . (|2.15p . in a more manageable form (see 
[M3| ). by introducing the change of coordinates 



JV-l , N-l 



Then 



a := s H kj, t — > r := t H kj. 

i=2 i=2 



*Txlt}=<* U\\boe*->) + 2 * 2 I dadK K)|V ?$0a 2 + D(K)+~X, (3.1) 
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$a [£] = (JV - 1) /_ d<xdrdK C(<r, K) £(r, K) (a 2 + r 2 + • r + £>(K) + A 

where 



K := k 



2, . . . , kjv-1 
iV-l 



|(a-,K) := i( cr 3_ V" ki,k 2 , . . . ,kjv-i ) , 

V m + 2 ^ ; 



D(K) :- 



m 



(m + l)(m + 2) 



JV-l \ 
(m + 3)£*f + 2£k i .k i J. 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



i=2 i<j 

Notice that -D(K) satisfies the bound (see (|3.49p in the following) 



m 



Yk$< d(k) < p^+E±A y k 2 

m + 1 ~ J 2 ~ (m + l)(m + 2) 4 



Last, setting 
and 

we obtain 



*:=^/D(K) + \p, r:=/D(K) + Aq, 

Q K (p) := (£>(K) + A) 3 / 4 £ fV^(K) + Ap, 



*^]=^K-a[|f[| 2 



L 2( R 3iV-3) 



f dK/D(K)+A Fi[Q K ] 

J M 3JV-6 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



where for any (>0we introduced the quadratic form in L 2 (M 3 ) 

^(F C ):=^(F 1 ) = |/GL 2 (M 3 ) jf d Pv ^+T |/(p)| 2 < ooj, F ( [f] := if iag [f] + ff [f] , (3.10) 



and 



F diag [;] = i dp 



m(m+2) 2 
(m+1) 2 



Ff [/] = " 1) / dpdq 



p 2 + c i/(p)r , 

/* (p)/(q) 
+ 1 2 + ^ttp • q + C 



(3-11) 
(3.12) 



Using the representation (13. 9p . (I3.10j) . (13. 11 H . (I3.12P we obtain the following estimate for $g, which 
is the crucial ingredient for the proof of Theorem 12.21 



Proposition 3.1 (Upper and lower bounds for $q) 
For any £ G £F(<1?^) we Ziaue 



l-A(m,iV)) < $m < (l + T(m,N)) 



s diag r 



where 



. , (iV-l)(m+l) 2 
r(m, iV) := = — arcsm 



y/m(m+2) 



m+1 



(3.13) 
(3.14) 
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The proof is based on a careful analysis of the form F^, £ G [0, 1], reduced to each subspace with 
fixed angular momentum I and it is postponed to the end of this section. First we introduce some 
useful notation and we prove some preliminary lemmas. 

For any / G L 2 (IR 3 ) we consider the expansion 



oo I 



/(p) = E E fim(p)Yr(0 P , <p p ) , (3.i5) 



1=0 m=-l 



where p = (p, P , 4> p ) in spherical coordinates and Y™ denotes the spherical harmonics of order I, m. 
We notice that / G S>(F\) is equivalent to 



oo I 



' — 1 - poo 

EE / dpp 2 ^/p^+i\f lm (p)\ 2 <oc. (3.16) 

1=0 m=-l ^° 

Moreover, we denote by P[ the Legendre polynomial of order 1 = 0,1,... explicitly given by 

Pl{y) = ¥v.ii {y2 - 1)l > ^[-M]- (3-17) 

In the first lemma we decompose in each subspace of fixed angular momentum I. 

Lemma 3.2 (Decomposition of F^). 
For f G @(Fi) we have 



oo I oo I 



F cif] = E E G <Aflm] =■■ E E (cf ag [flm] + G^l [flm]) , (3.18) 



1=0 m=—l 1=0 m=—l 



where for g G L 2 ((0,oo), p 2 \/p 2 + 1 dp) 



flf [g] := ^dp^J^^ + C \g(p)\ 2 , (3.19) 
GQ [g] := 2vr(iV - 1) fdp llq p 2 g* (p) q 2 g{q) fdy ™ — . (3.20) 

Jo Jo J -i p 2 + q 2 + T^nPiy + C 

Proof. For a given / G 3>{F{) we consider the expansion (|3.15p . From (|3.1ip we see that F c diag [/] = 

YlbLo Em=-I C diag [fim] ■ Concerning the off-diagonal term (|3.12p . we denote by 6 pq the angle between 
the vectors p and q and we consider the following expansion in Legendre polynomials: 

1 ^2Z + 1 f\ RM 



P 2 + q 2 + ^iPqcosOpq + C £^ 2 J_ x p 2 + q 2 + j^pqy + C 



oo „i 

E 2 ^ / d f 

7 n 1 



p ' (!,) Yl Yr%,<t> p )Yne q A q )- (3.2i) 



p 2 + q 2 + T^+rpqy + C 



In the last line we used the addition formula for spherical harmonics (see, e.g., |GR} Eq. (8.814)]). 
From we obtain Ff [f] = ZZo EL=-i Gfilflm]- □ 

In the next lemma we give a new representation of G^ which is particularly useful to control 
Gffi in terms of G^ for any ( > 0. 
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Lemma 3.3 (Estimates for Gf\). 
The form G^ can be written as 



G° ( S l [g] = Y, B ^ / dvv k e~^ / dpg(p) 



)p 2+k e-^ 2 



where 



Bik = ' 



otherwise . 



Moreover for any (>0 we have 



< Gfj [g] < GqH [g] fori even, 



(3.22) 



(3.23) 



(3.24) 



GZ[g]<G$[g}<0 for I odd. 



Proof. Using the expansion 



P 2 + q 2 + ^jpqy + C p 2 + q 



1 ^ / -2 



+ g 2 + C. 



and formula (|3.17p . we obtain 



k=0 



m+l)J J 



<W V, ' ^ \dyy k — l {y 2 -\) 1 



(p 2 + q 2 + (Y 



Integrating by parts I times we find 



cxj poo „2+k 



poo roo 

Gfi\9\ = Y, B hk® d P d q 
k=o Jo Jo 



Finally we use the identity 



A;! 



p_ g*{p)q '_ g(q) 

(p 2 +q 2 +C) k+1 ' 



(p 2 + q 2 + C) k+1 



dvv k e~^ + ^> 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



in (I3.28P and we obtain (I3.22|) . Let us fix I even. Then the integral in (|3.23p is different from zero 
only if k is even and this implies that is positive and the estimate (|3.24p holds. Analogously, 
when I is odd the integral in (|3.23p is different from zero only if k is odd and therefore G^J is negative 
and (^251) holds. " □ 



Now we study the form Go z = G + Gq ^ and we show that it can be diagonalized for each I. 

Lemma 3.4 (Diagonalization of Gqj). 

For any g £ L 2 ((0, oo), p 2 \Jp 2 + 1 dp) we have 



m + 1 



G$[g]= / dkS^ikf, 



(3.30) 
(3.31) 
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where 



g*(k) := -±= f dxe- lkx e 2x g(e x ) (3.32) 
V 2,-k Jm. 



and 

sinh f k arccos — %- 



/l sinn k arccos — ^rr 
dyPl(v)—- ^ _ J ■ (3.33) 

1 sin I Ftrc.c.nx — ^-r I Rmh(irk) 



sin ^arccos 7^75- J smh.(irk) 

Proof. The proof of (|3.30p is straightforward. To prove (|3.31|) we first introduce the new integration 
variables p = e xi and q = e X2 , so that the form reads 



G$[g) = 2*(N-1) [ dx 1 dx 2 e^g*{e x ^g{e^) C 

JR J-l 

= ir(N-l) [ d Xl dx 2 e 2xi g*(e Xi )e 2x2 g{e X2 ) I 

JR J-l 



p2x± i p2x2 _|_ 2y xi+xo 

m+1 (3.34) 

Pity) 

dy — TTf \~> — — ■ 

COSh(Xl — X2) + 



m+1 



The kernel in (|3.34p is a convolution kernel and therefore it can be diagonalized by means of the 
Fourier transform. Using the explicit Fourier transform of the kernel (see, e.g., [E]) we finally arrive 
at dS3U) . □ 

Owing to the previous lemma, the problem of finding bounds for the form Gq^ is reduced to 
finding bounds for the function Si(k). Taking into account the identity arccos z = § — arcsin z and 
the parity of Pi, we represent Si(k) as 



Si(k) 



-l sinh I k arcsin — ^~ 
-ir*(N-l) I dyP^y) ^ ' for I odd, 



i cos ( arcsin ) sinh (f A;) 



l cosh ( k arcsin 



v 



ir 2 (N -1)1 dyPi(y) ^ for { _ 



(3.35) 



cos (arcsin ^-j- J cosh k) 



It is evident from this representation that S[(k) is for any I > an even C°°-function of k with 
limfc-^oo Si(k) = 0. Before discussing upper and lower bound of Si(k), we shall prove the following 
elementary lemma. 

Lemma 3.5 (Taylor expansions of 5£ and 5|). 

For any fixed k > the following functions 

sinh f k arcsin — tt ) cosh f k arcsin — ) 
Sl(y) = -4 ^ > St(y) = (3-36) 



cos ^ arcsin J cos ^ arcsin 

have a Taylor expansion in the variable y £ [—1,1] with positive coefficients. 

Proof. Let V be the set of functions whose Taylor expansion has positive coefficients. First note that 
arcsin y, sinhy and coshy belong to V. The derivative is a linear automorphism of V and therefore 

d 1 1 m 
— arcsin y = — ^=^= = € V . 

dy J\ — yi cos arcsin y 

Moreover, V is invariant under dilations, multiplications and compositions of functions in V . Thus, 
S? and St belong to V. □ 
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In the next lemma we compute lower and upper bounds for Si(k). 



Lemma 3.6 (Bounds for Si(k)). 
For any k £l we have 



< Si(k) < 2vr 2 (iV-l) (m+1) arcsin 



1 



-47r(JV-l)(m+l) 



1 — \J m(m + 2) arcsin 



1 



m+1 



m + 1 
<S t (k) <0 



for I even (3.37) 
for I odd. (3.38) 



Proof. Let us prove (13. 38ft first. The upper bound follows from (|3.25j) and (|3.3ip . As for the lower 
bound, by means of (|3.36p we write 



Si(k) 



1T 2 (N-1) f + l 



sinh (ffc) 7_i 



(3.39) 



The first step is to prove that Si{k) is an increasing function of I for any fixed k. From (|3.39p . using 
(|3.17p and integrating by parts, we obtain 



Si+2{k) 



7T 2 (N- 1) 



sinh (|/c) 2 l+2 {l + 2)\ J_ x y dy 
ir 2 (N-l) 1 



^dy^(y 2 -^ l+2 s° k (y) 



sinh (f k) 2 l + 2 (l + 2)\ 

tt 2 (N - 1) 1 
sinh(ffe) 2 z + 2 (/ + 2)! 



l dy 



dy 



dy 



i 



(I + 2)(l + \){y 2 - l) l 4y 2 + 2(1 + 2)(y 2 - 



i QO 



d l S' 



dy 1 



(y) 



ir 2 (N-l) 1 
sinh (py ¥l\ J_ x 



^-D l %(y) + '- 



7T 2 (iV-l) 1+2(1 + 1) 



sinh(p) 2 i + 1 (/ + l)! J_ x 



dy (y 2 



ll CO 

1)1+1 _ D k 



dy 



f(y) 



Si(k) + 



ir 2 (N-l) 1+2(1 + 1) 
sinh(ffc) 2'+ 1 (! + l)! 7_i 



+i 



•/ co 



dy (y 



dy 



f(y)- 



(3.40) 



By Lemma 13.51 and taking into account that Z + 1 is even, we deduce that the last integral in (|3.40p 
is positive and therefore we conclude Si^ik) > Si(k). This means that it is sufficient to find the 
minimum of S x (k). From f)3.35|) we have 



S x (k) = -2tt 



\N-1) f 
Jo 



dy- 



y 



sinh ( k arcsin 



y 



m+1 



cos arcsin 



m+1 



sinh (f k) 



(3.41) 



We know that S x (0) < and lim / t_ s>00 S±(k) = 0. Moreover, the derivative of S x (k) does not vanish 
for k > 0, which follows from the fact that the derivative of the function 



sinh ak 



< a < b, 



sinhfr/c ' 

does not vanish for k > 0. Therefore, S x (k) is monotone increasing when k > and attains its 
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minimum at k = 0. Thus, 



Si(k) > St(k) > 5i(0) = -4tt(JV-1) / dyy 



arcsm 



m+l 



1 1 cos ( arcsin ^-j- 



-4vr(A - l)(m + l) 2 
-47r(JV- l)(m + l) 



dz z sin z 



(3.42) 



1 — y m(m + 2) 



arcsm 



1 



m+l 



and (|3.38|) is proved. The proof of (|3.37|) is completely analogous. In this case the lower bound 
follows from (|3,24p and f)3.31j) . Using the representation 



Si(k) 



7T 2 (A-1) f l 

sinh (ffc) 



(3.43) 



obtained from (|3.36p . one sees that Si + 2(k) ^ Si(k) and therefore it is enough to consider So(k). 
Since 



5 (0) = vr 2 (A-l) / dy — ^ = 2vr 2 (A-l)(m+l) arcsin (—^j) 

■'— 1 cos (arcsin —tt) \m + lj 

\ m+l I 



> 0. 



(3.44) 



lmifc-^oo So(k) = 0, and the derivative of So(k) does not vanish for k > 0, we deduce that So(0) is 
the maximum of So(k). □ 

Using the results of the previous lemmas we can finally prove Proposition 13. 1L 

Proof of Proposition HOI We prove first the estimate from below in (|3.13p . From Lemmas 
13.41 13.61 we have 



^i off [/] = E G °i ifim] + E G °Z l/u £ £gsjl/u = E / ^(ao|/L(*o 

7 ™ 7 7 ™ 7 ™ ■'K 



l,m 
I even 



l,m 
I odd 



l,m 
I odd 



l,m 
I odd 



> -4vr(A-l)(m+l) 



1 -Vm(m + 2) arcsin (— — -) E / dk \fL 



(3.45) 



-A(m 



where A(m,N) is defined in (|2.16p . Then, by ([22 



Ff ff [/] > -A(m,iV)E Gf^/zm] > -A(m,iV)E ag [/ hn ] = -A(m,iV) if ag [/] . (3.46) 



Lm 



Lm 



From the representation (|3.9p and from ()3.46p we have 

Site] > (1 - A(m, AT)) I dK \/ -D (K) + A F^ as [Qk] = (l - A(m, A)) S**^] 

JM.3N-6 



(3.47) 
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and the lower bound is proved. An analogous proof yields the upper bound in f|3. 13j) . We have 
Ff[f] < 2vr 2 (iV-l)(m + l)arcsin('-^ T W Jdk\f} m (k)\' 

l,m 

(iV-l)(m + l) 2 



y / m(m+2) 



■ arcsm 



(N-l)(m + l) 2 
< - — arcsin 



y/m(m+2) 

(N-l)(m + l) 2 
■sJm(m+2) 



^)E G o ias [/' 

l,m 

^)E G i ias [/' 



(3.48) 



rn 

diae 



which, together with (|3.9p . yields the upper bound for <£q. 



□ 



Let us briefly comment on the result of Proposition 13.11 By means of the elementary estimate 

(3.49) 



N-l N-l 



i=i 



i<j 



we find 



Kit] = aU\\h+m] < aU\\ L 2 + (l + T(m,N)) 

< / dk! • • • dk^Kk!, . . . , k*_!)| a «+2vr 2 (l+r(m,Ar))( 7 (m +;; j V fcf + A ) 



(3.50) 



and 



> / dki...dk JV _^(ki,...,k JV _ 1 ; 



i=l 



From (|3.50p . (|3.5ip . choosing A sufficiently large if a < 0, we conclude 

cfl-A(m,N))U\\ Hl/2 < <C||e|| ff i /a 



(3.51) 



(3.52) 



where c, C are two positive constants. Estimate (|3.52p implies that 1S finite for any £ £ 

and therefore our definition of the quadratic forms (|2.11j) and (|2.13p is well-posed. On the other hand, 
if in addition we assume A(m,N) < 1, from (|3.52p we conclude that <3?^ defines a norm equivalent 
to the if 1 / 2 -norm. This is the crucial ingredient for the proof of Theorem [ 

Proof of Theorem \2.2l From (j2.11j) and (I3.5ip we have 



^] > -A|Hli ? + A^[£] > -AiMii? 



+ Nj dki---dkjv_i |e(ki,... } kjv_i; 



a + 2vr 2 (l-A(m,A0) (J^^tf + A 

^ i=l 



N-l 



1/2 



(3.53) 
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for any ip G 3>(J- a ) and A > 0. Therefore, if a > the form T a is positive and if a < the lower 
bound (|2.20p holds. Let us now prove that T a is closed. We choose 



A > I 



if a > 0, 



4vr 4 1 



-A(m,iV)) 



if«<o, ( 3 - 54 ) 



and consider the form J„M := J r a ['i/ ; ] +A]|'0||?2 defined on QiiTa)- Let Upn} be a sequence in &(F a ) 



such that 



lim HV'n - VIIl? = , lim j£[V>n - Vm] = , (3.55) 

n— >oo i n,m— >oo 



where ip G if (M 3iV ). From the definition of #(.F a ) and J 7 ^ (see ([2+ll . (12+010 we have ^„ = 0*+<?a£„, 
with ^ G H{(R 3N ), e„ G F f 1/2 (M 3Ar - 3 ), and 

FabPn - 1pm] = ~ 0m] + ^*a[£n " U) • (3-56) 

This, together with (13.520 . implies that {cp^} is a Cauchy sequence in H^-(M. 3N ) and {£ n } is a Cauchy 

sequence in H^ 2 (M 31 ^ -3 ). Let us denote by <p> x and ^ the corresponding limits. From the explicit 
expression of the potential (|2.9p we notice that 

<ciien i? , (3.57) 

where c > 0. Hence, 

lim [|^-(0 A + 0aOIIl?= lim ||(^-^ A ) + g A (^-e)llL f2 =0. (3.58) 

Since the limit of the ip n 's is unique, xp = (p x +Q\£. Therefore xp G 9{F a ) and lim^oo F x [i)-ip n ] = 0. 
This shows that the form J^ x is closed and a fortiori J 7 ^ is. □ 



4 Unboundedness from Below of T a 

This section is devoted to the proof of Theorem 12.31 As we shall see, what makes an instability 
condition hard to prove is the restriction to antisymmetric wave functions. 

In fact, the proof relies on the explicit evaluation of the charge form <I>^ on a trial function, i.e., 
a convenient sequence of charges with energy going to — oo. Identifying one such sequence is easy 
when N = 2 because 3> x 18 m practice the same as the reduced form F% - see (14. 3|) below - and the 
analysis performed in Section [3] suggests that a convenient Q n (p) has to be chosen in the subspace 
with angular momentum I = 1 and such that in the position representation it becomes peaked at 
the origin as n — > oo (i.e., two identical fermions coming arbitrarily close). 

When N > 2, on the other hand, a natural trial function satisfying the antisymmetry constraint 
would be the Slater determinant of N one-particle charges, one of which is Q n itself. A convenient 
choice is driven by the physical idea of a A-body configuration that contains precisely the (2 + 1)- 
body structure minimizing the energy with A = 2, whereas all remaining particles are placed far 
away in space so that there is no or negligible interference with the two-body state. This results in 
a A-particle Slater determinant between Q n and A — 1 copy of a different component (see (|4.19p 
below). The fermionic character of the trial function is thus fulfilled by construction and optimising 
the choice of the second component produces only higher order symmetry correlations. 

Throughout this section we assume that A is a positive number such that G\ < A -1 < C2 for two 
finite constants C\,C-i < 00, which in particular will allow us to incorporate error factors proportional 
to A -1 into a constant C. 
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Proof of Theorem \2.3[ In order to prove instability of the form T a , it is enough to produce a sequence 
of normalised charges £ n G Hj' 2 (WL 3 ( N ~ 1 )) such that 



lim $^[£ n ] = -co, 

n— >oo 

since the sequence of states GxCn then satisfies 

Fa[G\(,n] = -A ||£/A£n|li2( K 3iV) + .V«i>;Vi„^ 



->■ — oo. 



(4.1) 



(4.2) 



Case N = 2 . The result was already proved in [FTj . but we repeat here the argument for we use a 
slightly different trial function that turns out to be useful in the general case N > 2. Owing to (|3.9p . 

^K]-a[|eili a( BS ) = V^Fi[Q] (4.3) 

where Q(p) = A 3 / 4 £(\/Ap) (recall (|3.8p ). Then we only need to produce Q n G L 2 (M 3 ) such that 
lirrin^oo Fi[Q n ] = — oo. Note that no constraint is imposed on the symmetry properties of Q n . In 
fact, according to the discussion of Section [3] (see (|3.30j) . (|3.3ip and (|3.42p ). we can take each Q n in 
the subspace with angular momentum I = 1 and such that the support of its ft— transform defined in 
(|3.32p gets concentrated at the origin. Explicitly, we choose 



Q n , 7 (k) := n- 3 / 2 Q 7 (n- 1 fc)y 1 °(i? fc ) ) 

Q 7 (p) := 7T- 1 /4 &r7 V2 p -l exp |__^ j exp |_1 7 2 (lQgp) 2 j @(p _ 1} 



(4.4) 
(4.5) 



where 7 G (0, 1) is a variational parameter, is the Heaviside function, i.e., O(p) = 1 if p > and 
otherwise, and c 7 is a normalisation constant. By direct computation, 



IIQn. 



7llL 2 



dpp 2 \Q y (p)y 



A r dt exp {-t 2 + § - ^} = |c 2 (l + erf {^}) 
•/ 

Imposing Q nn to be normalised yields 

l<c 2 = 2 [l + erf{^}]~ 1 < l + CTexpj-^}. 
(see, e.g., |AS} Eq. (7.1.13)]). It is also useful to compute the following integral for a > — 1: 

J dpp 2+a |g 7 (p)| 2 = ^ 7 c 2 exp{- 3 i 7 }y o dt exp{- 7 2 t 2 + (l + a)i} 

= ±c 2 (l + erf { ^}) exp {^(2a + a 2 )} < (1 + C 7 ) exp {^(2a + a 2 )} , (4 

where we used 



(4.6) 



(4.7) 



ic 2 ( 1 + erf 



1 + erf {^) 2 



2l 



l + a 

27 



df 



vr /j_ 

2 7 



e"* 2 < 1 + C7- 1 exp {-^} < 1 + 0( 7 ). 



Using the decomposition (|3.18p . as well as the scaling law of Q n ^ with n, we have 



-^1 [Qnr, 



G 



diag 



n- 3 / 2 Q 7 (n- 1 p)l + Gf x In-^Q^p)] = n \g*$ [Q 7 ] + G° ff 2j [Q 7 ]l . (4.9) 
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We estimate the diagonal term in (|4.9p as 

a^ m < 2 ^V!^) (1+o( „-i )) || »||^ 



(4.10) 



£2T2 v^i exp {^} [1 + 0(0 + 0(7)]> 

where we used 

||Q 7 fc = / dfce 4fc |Q 7 (e fc )| 2 = /° dpp 3 |Q 7 (p)| 2 <(l + C 7 )exp{ A ). (4.11) 

As for the off-diagonal term in (|4.9p . 

G^ a>1 [Q 7 ] = Gg[Q 7 ]+7J, (4.12) 

where 



\Tl\<C ds / dt 



s 2 + t 2 + -^s • t + n' 2 ) 1 - fs 2 + t 2 + -^s • t N 

m+1 / \ m+1 



|Q 7 ( S )||Q 7 (i)| 



<Cn 2 I ds [ dt i%MJ%ffl < Cn -2 / ds |q (s) |2 r dt r 2 <0(n- 2 ). (4.13) 

Moreover, 

Gfi[Q 7 ] = f dk S^Q^k)] 2 
Jr 



< S 1 (0)\\qU\%+ / dk (5i(fc)-5i(0))|Q» 



< 51(0) exp {^} (1 + 0( 7 )) + C^dA: v^|Q 7 (&)| 



(4.14) 



2 " 2 ^TT^ A(m ' 2) exp { ^} (1 + ° (7)) + c I dk vWt(*)I 



where we used (|3.42[) and the elementary estimate Si(k) — S\(0) < Cy|fc|. To estimate the last 
integral in (|4.14|) we observe that 



Q^(jfe) = Tr-^^^expj-g^} -j= J dx expj-ifoc- \x 2 + xj 



7r- 1 / 4 c 77 - 1 / 2 



exp {-^} (exp {-£ - ^ + ^} - ^(1 + r(*))) , 



(4.15) 



where 

\ — ik * / \ i 7 2 , 

Z = W |r(z)|£ TO- (4 - 16) 

Therefore, 

^|Q 7 (,)| 2 < ^expj^^^exp^j^exp^}^^! 



< CV7exp{^}(l + V7exp{-^r}). 
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Using pO| . flU . (jTOD . gUD, (HTTjl in flM]) we finally obtain 



Fi[Q n , 7 ] < 2 ^ m{m J 2) ne Xp \A\ [l - A(m, 2) + + O^ 1 )] ► -oo, (4.18) 

if A(m, 2) > 1 and 7 is taken small enough (independent of n). 

Case N > 2 . As mentioned at the beginning of this section, this case is more complicated, for the 
trial sequence £ n must be antisymmetric under the exchange of any variable, i.e., £ n G Lf(M. 3N ~ s ), 
and at the same time we want £ n (ki, . . . , k/v-i) to behave like Q n ,7(ki) once the other degrees of 
freedom are traced out. Looking for £ n matching these two requirements is an example of the well- 
known represent ability problem (see, e.g., [LSJ), i.e., the search for sufficient conditions to impose on 
a one-particle density matrix so that it can be obtained as the reduced density matrix of a fermionic 
many-body state. We remark that the solution is known only in some special cases and is non-trivial. 
Our choice here is a trial state that is as close as possible to an uncorrelated state, which is given by 
an antisymmetric wave function containing Q nrr Explicitly, 



£n(kl 



Qn,-y(kl) 



1/3,2 (kl) 
l/3,2(k 2 ) 



'/3,JV-l(kl) 
'/3,iV-l(k2) 



Qn,7 (k 



S /3,2(k 



N-l, 



^/3,iV-l(kjV-l 



(4.19) 



where Q n ^ is defined in (|4.4p . < (3 <C 1 is another variational parameter, 

S w (k) := (4vr)- 1 / 2 /3- 3 / 2 -(/T^) exp{il<p k } , 
I G N, H G C^°(]R + ) is real- valued, with support in (0, 1), and such that 



dk k 1 



\k) 



1. 



(4.20) 



(4.21) 



Note that, since the two functions Q n , 7 and Z > 0, are orthonormal by construction, the 

function (14. 19ft belongs to L 2 (M 3 ^ -1 ^) and is normalised. Moreover, the supports of Q 7 and S do 
not intersect, which implies that the supports of Q n , 7 and j are disjoint as well, provided /3 < n, 
which follows from the assumptions on f3. 

We can now evaluate We start by estimating the diagonal part. Using the exchange 

symmetry and the definition of L\ in (12.80 ). 



a + 



1 



(N - 2)! 7 R 3(jv-i) 



dkxdK L A (ki,.. 



, kjv-i 

N-l 



^ [] sgn( ff ) S gn(r)3^,(k (7(0 )H Ai (k r(i) ), (4.22) 



a,TEV N -i l,j=2 



where Vn-i is the group of permutations of N — 2 elements 2, . . . , N — 1 and sgn(cr) denotes the sign 
of any a G Vn- All the other terms vanish because of the integral of the product Q nj7 (kj)Hg(kj), 
which is pointwise zero thanks to the disjoint supports of the functions. Extracting the main factor 



m(m+2) ,2 1 \ 
(m+l)^ K l + A ' 



and bounding the rest by means of the inequality 

Va + b < y/\a\ + <s/\b\, 



for a + b > , 



(4.23) 
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we obtain 



N-l 



m(m+: 
(m+1) 



2m 



-1/2 



m(m+2) \p,2, 

(m+1) 2 * (m+1) 2 



i=2 



i>i 



i<i<i 



1/2 



(4.24) 



N-l 



< + All + C N £ fa + Vhk?j } ■ 

The diagonal term can be estimated as 

- a < 2vr 2 / dki • • • dkjv-i J^^Tl |Qn, 7 (ki)| 2 n |S/j,o(ki)| S 

N-l 

IQn, 7 (ki)| 2 H |S AO (kO| 2 , (4.25) 



7V-1 



+ Cat / dki • • • dk^v-i 

/ R 3(iV-l) 



L A ( kl ,K)-2.yi^^ + A 



Z=2 



thanks to the orthogonality of functions H^j and E^/ for I ^ V . Thus, by (|4.24p . ^ Q A g [£n] — a is 
bounded from above by 



2vr 2 / dk 1 v/^^fc? + A|Q n , 7 (k 1 )| 2 / dK 

J R3 V (rn+i, J R 3(N-2) 



l + C N [k 2 + Vhk, 



N-l 



II l H Ao(k/)P 



1=2 



< 27r 2 n V /T ^ 1 +2) ( X + 0(„-i)) J dkt k\ [l + Cjy (/? + \A^T)] IQ-yC^i 

< 2vr 2 n^^ exp {^} [l + C N (^exp {^} + 7 + fi + 



where we used (|4.8p . We now compute the off-diagonal term (recall (|2.15|) ). Owing to the exchange 
symmetry, the pre-factor N — 1 cancels with the normalisation factor of £ n and 



®f[U = [ dsdtdKG A (s,t,K)(Q* (s)Q n , 7 (t) 

JR3N { 



(N 



1 7V_1 1 

-— y n sM°)sm(T)zh( k °(i)) E eAKu)) +^ (4-27) 



<T,re7>jv_i ZJ=2 



where 7?. contains some remainder terms. We estimate the leading term (first term on the r.h.s. of 
(|4.27p ) from above by 



i> i dsdtG A/ „ 2 (s,t)Q 7 ( s )Q 7 (t)y 1 ( 1 ? s )y 1 (^) 

+ C N [ dsdtdKG A (s,t,K)G A (s,t) [kl + {s + t)k 2 ] |Qn, 7 (s)| |Q ftl7 (t)| TT l^k,)) 2 , (4.28) 



AT-1 



Z=2 
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where in the first term we replaced G\(s, t, k 2 , . . . , k-N-i) with G\(s, t, 0, . . . , 0) =: G\(s, t) and we 
exploited the orthogonality of functions Eri and Epp f° r I ^ V • The first term in the expression 
above was bounded in (|4.14|) . Using 



G A (s,t,k 2 , . . . ,kjv_i) < 



m 



m + 



-{s 2 + t 2 )+\ 



(4.29) 



and the elementary inequality x + y < V x 2 + \\J y 2 + 1, x,y > 0, the second term in (|4.28[) can be 
estimated as 



N-l 



dsdtdK k2 + (S ± f 2 |Q n , 7 (s)| |Q B , 7 (t)| [J |S/j,o(k, 
[s 2 + t 2 + A^] 



< C^n -1 / dsdt 



/3 2 + n/3(s + t) 
[ s 2 + t 2 + \m±L n -2 



< C N n- 1 (0 2 + n/3) 



1=2 

2 |Q 7 (s)| |Q 7 (t) 
1 



(4.30) 



sup / ds- . 



< 



[s 2 + t 2 + X 
C7v(/3 2 + n/3)(l+ 7 )exp{^} . 



dp (p 2 + l)|Q 7 (p) 



The rest 1Z in (|4.27p contains several terms but it is not difficult to see that most of them vanish 
because of the disjoint supports of Q Ujl and Ep. What remains is 



[(N - 3)!]- 1 / dsdtdk 2 • • • dk N -i G x (s, t, k 2 , . . . , kjv-i) El ^Ep^t) |Q„ l7 (k 

JR3N 



2j| 2 X 



N-l 

x II s g n ( CT ) s g n ( r )^/ k <x(o) H /3.i( k ^(i)) 

a,TeVN-2 i,j=3 

3 



(4.31) 



< Cjv / dsdt |H^,i(s)| |H^,i(t)| < Cn/3*, 



where we exploited the exchange symmetry again, as well as the properties of H (in particular 
supp(S) C (0,1)) and (|4.29 jl . It is understood that the sum over permutations as well as the 
following factor in (|4.31j) is not present when N = 3. Putting together (14.271) . (j4.28j) . (|4.30h and 
(|4,3ip . we obtain 

*f [£n] < nGf /n2jl [Q 7 ] + C N (f3 2 + n/3) exp {^} , (4.32) 

and finally 



^]<2vr 2 nV^±^exp{ i7 .}{l-AK2) 

» _1 + V7 + Vrf^{jjfr} + (P~ 1 P 2 +P)exp{$*}] }• ( 4 - 33 ) 



+ C N 

By assumption, 1 — A(m, 2) < and we choose (3 <C as, say, 

/3 = n - 2 . Thus, we can always find 

some small 7 = 0(1) > 0, such that 

*£Kn] ► -00 (4.34) 

n— >oo 

which concludes the proof. □ 
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Appendix 

Here we describe the formal procedure for the construction of the quadratic form T a - We start from 
the Hamiltonian fj2.4|) written in the Fourier space 



N 2 \ 

m)(k 1 ,...,k N ) = [ ^fc l 2 +^-^k i -k,)v(k 1) ...,k iV 

i=l m+ i<j ' 

+ -^- / ds^(ki,...,k t-i)S,ki + i,...,kjv) (A.l) 
{2-kY J r3 



and we consider the corresponding quadratic form, regularized by means of an ultra-violet cut-off 

N 



FM-= I dk 1 ...dk A r{f^fe i ? + r ^— J>-kA i>{k u ...,k N ) 

+ 4t^tE / dk 1 ...dk NXR (k i )r(k 1 ,...,k N )x 

( 27T ) ~1 Jr 3n 



x / dsxii(s)^(ki ) ...,k i _i ) s,k i+1 ,...,kjv). (A.2) 



Here ip E i7 1 (M 3Ar ), X-r( s ) is the characteristic function of the three-dimensional ball s < R, and a 
is a parameter that has the role of a renormalised coupling constant. Note also that we introduced 
in fi a dependence on R: the choice of such an explicit dependence will be the main content of the 
renormalisation procedure. 

We now define the "surface charges" £f £ L 2 (R 3N ~ 3 ) associated with tp G H 1 (R 3N ) as 

If (ki, . . ■ ,kjv-i) := ^"'^P / dsxii(s)^(k 1 ,...,k i _i,s,k i ,...,kjv_i), (A.3) 

and the corresponding "volume charges" pf (ki, . . . , kjy) := Xn(ki) £f (K*)- Further, we introduce 
the "potential" 

TV 

S^Oci,...,^) :=^G A (k 1 ,...,k JV ) Xi? (A; i )ef(K 4 ), (A.4) 
i=l 

where Ga is defined in (|2.7p for any A > 0. Setting 

0« := ^ - Q^R (A.5) 

we have 

= [</>f] + A ||" 2(R3N) - A ll^ll^^) + A , (A.6) 
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with Fo[<j)\ := Hq\4>), and 

JV 



i=l 



^' A [?]:=-E / dk 1 ---d^ NX R(k i )i*(K i ) i>(k 1 ,...,k N ) + G x (k 1 ,...,k N )£ i (K i ) 



V / dki- • • dkjy xdh) i*(ki)G x (ki, ...,k N ) XR (kj) ijiKj) . (A.7) 



In the limit R — > oo we assume that pf~, ^f 1 — > £j. Moreover, we extract from the diagonal part of 



(|A.7p only the terms not vanishing in that limit 

(2vr)3 



TV 

E 



dK ?: 



R3iV-3 



£i(Ki) 



/x(a, i?) 

/V 



dki XR(^)G A (k! 



dK, : 



&(K;) 



(2vr) 3 



4vri? 



+ 2vH 



m(m + 2) 
(m + 1) 



V/c 2 I 2m Vk k I A 



1/2 



+ 0(1) 



In order to remove the cut-off one is thus forced to set fi — > as R — > oo and, although several 
choices are allowed, we set 

(2vr) 3 



jj,{a, R) :- 



(A. 



AttR + a 

this way canceling the singular term proportional to —4irR contained in the expression above. 

We can now remove the cut-off taking the limit R — > oo and so recovering the expression (|2.1ip . 
Note that we exploit at this stage the fermionic symmetry, which in particular implies that all charges 
can be expressed in terms of a single function £, i.e., 



&(xi, . . . ,xat_i) = (-l) l+1 £(xi,...,x/v_i), 



(AS 



and £ itself is totally antisymmetric under exchange of coordinates. This in turns implies that the 
sign in front of the off-diagonal term is the opposite than in the bosonic case, implying a completely 
different behavior of the ground state. 
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